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A Comment on the Derivation of the 
Langmuir Isotherm for Adsorption 

Kinetics* 

C.L. Frenzen 

Department of Mathematics 
Naval Postgraduate School 
Monterey, CA 93943, USA 



Abstract 

The justification for the use of the Langmuir isotherm, basically 
an equilibrium relation, to relate the subsurface concentration C(t) 
at time t to the adsorbed surface concentration T(t) depends, in many 
adsorption kinetics models, on the presence of two very different time 
scales. Though this was implicit in the recent derivation in [1], it was 
not specifically pointed out. The second dong’ time scale is partic- 
ularly important because it governs the time evolution of C(t), the 
subsurface concentration, and helps to determine when the Langmuir 
isotherm approximation is valid. Here the Langmuir isotherm is de- 
rived mathematically from the adsorption kinetics equations governing 
the surface concentration, the subsurface concentration and the bulk 
phase adjacent to the surface. 



The dependence of T 00 , the equilibrium concentration of a surface-active 
component at an interface, and C°°, the corresponding bulk phase concen- 
tration, are often related by the Langmuir isotherm (Langmuir [2]): 

* International Classification for Physics numbers 5100, 5110 
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(1) 



poo 7rc~ 

1 +7 C°° 

In ( 1 ) f represents the saturated surface concentration while 7 is a constant. 
Frequently the assumption is made that at any time t the adsorbed surface 
concentration T(t) and the ‘subsurface’ concentration C(t) (ie., the concen- 
tration in a thin layer of bulk phase adjacent to the surface) are related by 
the equation 



m = 



7f C(t) , 

1 + 7 c(ty 



(2) 



see the references given in [ 1 ]. Equation (2) is often used as a boundary 
condition to determine the surface concentration from the subsurface con- 
centration at time t. Some of the the inadequacies of making the intuitive 
assumption in ( 2 ) and the desirability of replacing it with a more fundamen- 
tal description of the Langmuir isotherm have been discussed in [ 1 ]. We also 
note that in many time dependent problems it is important to establish how 
initial values of the surface and subsurface concentrations evolve over time to 
eventually satisfy (2). This time evolution can greatly affect surface tension. 
For example ( see [3], p. 62 ), at the surface of ordinary tap water the surface 
tension immediately after formation of the surface is close to the value for 
pure water, but usually rapidly falls to something like half this value due to 
adsorbed contaminants. 

The key to the validity of ( 2 ) lies in the presence of two very different 
time scales: £ p , the characteristic time defined by the rate constants for 
adsorption and desorption of subsurface surfactant into the surface, and t s , 
the characteristic time defined by the rate constants for passage of subsurface 
surfactant into and out of the bulk phase. The letters p and s stand for pre- 
steady state and steady state respectively. When t p <C t s the relaxation in 
surface concentration for a given change in subsurface concentration is very 
rapid. The time scale t s is particularly important because it controls the 
time evolution of C(t), something which could not be determined from the 
derivation given in [ 1 ]. The purpose of this letter is to extend the model given 
in [ 1 ] to take into account the passage of subsurface surfactant between the 
subsurface layer and the bulk phase, thereby deriving the Langmuir isotherm 
and conditions for its validity from an adsorption kinetics point of view. 
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As in [1], let J + be the rate at which subsurface surfactant is adsorbed 
and J~ be the rate at which adsorbed surfactant is desorbed. Assuming no 

evaporation, we have 

d 

— T = J + -J~. 

dt 

We assume, as in [1], that 

J+ = kiC(t) fl - , J~ = k 2 V(t) 

where ki and k 2 are rate constants. Consequently 

d 



dt 



T(t) = k x C{t) 



(■-¥)- 



k 2 r(t). 



( 3 ) 



Let I< + be the rate at which bulk phase surfactant passes into the subsurface, 
and K~ the rate at which subsurface surfactant passes back into the bulk 
phase. Then 

4-C(t) = J~-J+ + K+ - I<~. 

dt 

We will assume that if C°° is the equilibrium bulk phase concentration, then 

K + ~ I<~ = k 3 (C 00 - C(t)) 

for some rate constant & 3 . Therefore 

jC(t) = (l-Wt'j- k 2 r(t)) + t 3 (C“ - C(t)). (4) 

It is equation (4), governing the time evolution of C(t), which was not present 
in [1]. As initial conditions to (3) and (4) we assume 

r(o) = f 

C(0) = C°°. (5) 

In many cases f may be zero, but it need not be. 

To understand more fully (3), (4) and (5) we shall now nondimensionalize 
them by letting 

V(t) = C°° T*(C) 

C(t) = C°°C*(C) (6) 

( = f 1, 

k 3 
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where starred quantities are dimensionless. We have chosen l/k 3 , the char- 
acteristic time for passage of subsurface surfactant into and out of the bulk 
phase, to nondimensionalize t. Thus t = t a t m where t 3 = 1/&3. As in [1], we 
shall let p = fci /&2 and for purposes of discussion here we shall assume that 
p is an order one quantity. The other cases discussed in [1] can be easily 
adapted to this one. 

By our assumption on p, 1 / is a, characteristic time for adsorption and 
desorption into the surface, and we therefore let t p = 1/&2- Defining 



£ = 




h 

V 



we shall assume that 

£ < 1 . 

Physically this corresponds to assuming that the characteristic time scale 
t p over which changes between the surface and subsurface concentrations 
occur is much shorter than the time scale t s over which changes between the 
subsurface and bulk phase concentrations occur. 

Inserting (6) into (3), (4) and (5) and dropping the stars for convenience, 
we find their dimensionless form to be 

€j t C(t) = - (pCft) (i - Slr(o) - no) + 1 (! - C(0) (8) 



r <°> = ^ 

C(0) = 1. (9) 

Note that all dependent and independent variables in (7), (8) and (9) are di- 
mensionless 0(1) quantities. The system (7), (8), (9) is a singular perturbation 
problem, characterized by the presence in (7) and (8) of a small parameter 
multiplying the highest derivatives. We cannot let e — » 0 and still hope to 
satisfy the initial conditions in (9). The resolution of this problem is found 
in the existence of a ‘boundary layer’ in time where T(<) and C(t) change 
so rapidly that one cannot neglect the derivative terms in (7) and (8). This 
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region corresponds, as we shall see, to the rapid relaxation of the surface and 
subsurface concentrations from their initial values. After this rapid relax- 
ation, or pre-steady state period, the surface concentration T(t) is in steady 
state with respect to the instantaneous value of the subsurface concentration 
C(t). This yields the Langmuir isotherm given in (2). 

Away from the region of rapid transition, ie., in the steady state, we look 
for the so-called ‘outer’ solutions in the form of asymptotic expansions in e : 

ro^) = ro (*) + e r, (* ) + . . . , c outeT {t) = c 0 {t ) + £ c, (<) + ... . (io) 

Substituting (10) into (7) and (8) then yields to 0(£°): 

pCo(t) ^1 — ro(<) = 0, 

or 

r o (0 = . — . 

r + pc~c 0 (t) 

This is, of course, the dimensionless Langmuir isotherm given in (2) with 
7 = p/T . The 0( e°) surface concentration r 0 (t) is in steady state with 
the instantaneous value of the 0(e°) subsurface concentration Co(t). However 
(11) yields no information about the time dependence of Co(0> and we cannot 
possibly satisfy the initial conditions in (9). To discover the former, we add 
(7) and (8) together and find, after cancelling a common factor of e, that to 

0(e°) , 

J t (r o (0 + c„(0) = 1 - C 0 (t). 

This coupled with (11) then leads to the following single differential equation 
for C 0 (t): 

Jt c «) “ • < 12 > 

1 T 7 7T 

(r + pC°°c 0 {tj) 

Since C(t) is a concentration, we will always require Co(0) > 0. We 
should not necessarily require C o (0) = 1, however. We are discussing the 
outer solution, and rapid changes may already have taken place in T(^) and 
C(t) making them differ from their initial values. We will shortly discover 
how to choose Co(0) by ‘matching’ it with the rapidly varying or ‘inner’ 
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boundary layer solution. C o (0 = 1 is a solution to (12) and it is not hard 
to see that when Co(0) > 0 all solutions of (12) tend to one as t —* oo. For 
example, if 0 < Co(0) < 1, (12) implies that 



so that 



0 < 



j, c ^ 

1 — C 0 (t) 



<1, 



Co(0) < C 0 (t) < 1 - c-‘(l - Co(0)) 
for all t > 0. Using (12) and (13) we then find that 

d _ . „ 

7t c ° {t) - 

‘ T (f + />C“C o (0)) 2 



1 ~ Co{t) 

Pp 



and integrating this inequality gives 



1 - C 0 (t) < (1 - C o (0))e- Alt , 

where 

Ai = hr . 



Combining (13) and (14) together then gives 

(1 - Co(0))c"‘ < 1 - C 0 (t) < (1 - C o (0))e- Alt 
if 0 < Co(0) < 1 . In the same way, if Co(0) > 1 we find that 
(C o (0) - l)c-‘ < C 0 (t) - 1 < (C o (0) - l)e" A21 



1 

ry 

(f + pC">C 0 (oj) 



where 



Ao — 



1 + — 



r V 



(f + pC-)' 



(13) 



(14) 



(15) 

(16) 
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We note that (12) can be solved implicitly by separation of variables, but we 
find that (15) and (16) supply more information. Equations (15) and (16) 
imply that all solutions of (12) with Co(0) > 0 tend to one as t — > oo, at a 
rate controlled by the characteristic time t 3 = l/k 3 . In dimensional notation 
recall that this means that the 0(e°) outer solution C°°Co(t) tends to C°° 
as t — ► oo, and hence from (11) that the 0(e°) outer solution C°°ro(2) ( in 
dimensional notation ) tends to 



fpC 00 

t + pC°° 



as t — > oo. This is precisely the form of the Langmuir isotherm given in (1) 
with 7 set equal to p/T. 



clT* dC 

What about the ‘inner’ or boundary layer solution in which — and — — 

dt at 

change rapidly? To see this rapid change we introduce a ‘stretched’ dimen- 
sionless time t defined by t = -. Since —z. = e — , when derivatives with 

e dtdt 



respect to t are 0(1), derivatives with respect to t are 0 and rapid 



changes are occurring. Note that t = k 2 td = td/tp where id denotes dimen- 
sional time. Hence t measures dimensionless time at a rate characteristic of 
the rapid processes of adsorption and desorption into the surface. 

Using t instead of t in (7) then gives 



jT(t) = pC{i) (l - &(?)) - r(i) (17) 

fx(i) = - (pC(t) (l - ST(()) - r(i)) + £(1 - C(t)). ( 18 ) 

These are the differential equations for the ‘inner’ or boundary layer solution 
describing rapid changes in V and C. Again we look for solutions in the form 
of asymptotic series in e: 



r er (f) = f 0 (i) + £f 1 (i)f ... 

C inner (i) = C 0 {t) + eC 1 (t) + .... 
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Substituting these into (17) and (18) and using the initial conditions in (9), 
we find the 0(e°) ‘inner’ problem to be 

d / c°° \ 

„(!) = pCo(i) (1 - — fo(f)J - r o(i) (19) 

jCoit) = - (pC„(i) (l - yTo(<)) - Mi)) (20) 

MO) = ^j,Co(0)=l. (21) 

Adding (19) to (20) and applying (21), we conclude that 

f 0 (f) + C 0 (i) = B (22) 

where the constant B is defined by 

B = 1 + (23) 



Thus, on the boundary layer time scale t, the sum of the surface and sub- 
surface concentrations is conserved. Using (22) to eliminate Co(t ) from (19) 
then gives 

im = p (M?) - f ,) (Mi) - f 2 ) (24) 



where 



I\ = 



C 



F (i + -) + B + 



C 



r „ U \ 4 BY 

(l + -) + B) - 



C c 



r, = 



c 



r i 

(i + -J + B- 



' r „ U n \ 4BT 
^(l + -) + i?) 



(25) 



Both fi and f 2 are real and positive, and since 

\ 2 



T „ 1. \ 4BT 

— (l + -) + 5 



C 
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we find from (25) that 



f 1 

— (! + -) + B + 

C°° v p J 




2 




f , 1 , 


r , i N 


c ~ (1 + ; ) + B - 


c~< 1 + p)- B 



It follows from (26) that 

fi > B 
f 2 < B. 



( 26 ) 



(27) 



Note that fj and f 2 are equilibrium solutions of (24). If f o (0) ^ fi and 
f o (0) 7 ^ r 2 , (24) may be solved using the initial condition 



f o (0) = 5-1 



from (21) and (23) to yield 



fo(f) = 



r, - r;Z 0 e»‘ 

1 - Z 0 e°‘ 



(28) 



where a > 0 is given by 



a = p- 



C° 



f \ 




and 



Z 0 = 



f o(0) - f ! _ B - 1 - f 1 

B - 1 - f 2 



(29) 



r„( 0 ) - r 2 

From (28) and (29), r o (t) blows up in finite time if and only if 0 < Z Q < 1 , 
or if and only if < B — 1 . By (27), however, this can never happen. 

It follows from (28) and the above remarks that, whatever fo(0) is, all 
solutions f 0 (t) — * r 2 as i — » oo. Hence, from ( 22 ), Co(t) —* B — T 2 as t 
Note that as £ — > 0, t — ► oo when t differs from 0 by an 0(1) amount. 



oo. 
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We shall now ‘match’ the inner and outer solutions together for both the 
0(e°) surface and subsurface concentrations by requiring that the limit of 
the inner solution, as t — * oo, ‘matches’ or is equal to the limit of the outer 
solution as t — * 0. This implies that 

limT 0 (t) = lim f 0 (<), lim C 0 (t) = lim C 0 (t) 

<— 0 t— oo <— 0 t— CO 

or , using (11) and (28) 



f>Co(0) 

f + P C°°Co(0) 



— f'2,C'o(0) — B — f 2 . 



(30) 



From (27), Co(0) > 0, and we see that the matching conditions in (30) 
determine initial conditions for the outer solutions r 0 (t),Co(t). 

The final step is to form the ‘composite’ or uniformly valid solution by 
adding the 0(e°) inner and outer solutions together and subtracting the 
common parts as given in (30). This yields 

r(f) = r o (0 + fo(i)-f 2 + O(e),C(0 = C 0 (t) + C 0 (i)-(B-t 2 ) + 0(e). (31) 
Recalling that i = j, we can write (31) as 

r (t) = r o (0 + f 0 —t 2 + O (e) , C (t) = Co(t) + Co — (B — f 2 ) + 0(e). 

( 32 ) 

From (32) we see that the solutions for T(t) and C(t) each consist of an outer 
solution plus a boundary layer ‘correction’. Equation (32) gives the uniformly 
valid composite solution to 0(e) for the dimensionless system given by (7), 
(8) and (9). 

Initially, by (21) and (30), Fo(0) = f 2 and fo(0) = B — 1 = T /C°°. Then 
as t increases Io(t) rapidly approaches F 2 .This domain of rapid change is 
the pre-steady state. After the pre-steady state the steady state begins. In 
this regime To(0 ~ T 2 and, to 0(e), T(t) = T 0 (t), where the outer solution 
T o (0 is given by (11), with C 0 (t ) satisfying (12). By (15) and (16) the outer 
(dimensionless) solution C 0 {t) tends to one for long times at a rate controlled 
by < s .This means that T(t) tends ( in dimensional form) to the expression in 
(1) with 7 = p/T. We note that by using results in [4] we can prove that the 
composite solution in (30) is a uniformly valid solution to 0(e) of the system 
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given in (7), (8) and(9). The reader can find further information on solution 
methods for problems of this type in [5] and [6]. 

To conclude, the validity of the derivation of the Langmuir isotherm given 
here depends on the presence of two different time scales t p and t a , with 
t p /t 3 = £ < 1. These time scales characterize the pre-steady and steady 
state time scales respectively. It is after the relatively short pre-steady state 
period that the 0(£°) surface concentration Fo(t) is in steady state with the 
instantaneous value of Co(t), the 0(e°) surface concentration, yielding the 
(dimensionless) Langmuir isotherm given in (11). 
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